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Cusp-scaling behavior in fractal dimension of chaotic scattering 
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A topological bifurcation in chaotic scattering is characterized by a sudden change in the topology of the 
infinite set of unstable periodic orbits embedded in the underlying chaotic invariant set. We uncover a scaling 
law for the fractal dimension of the chaotic set for such a bifurcation. Our analysis and numerical computa- 
tions in both two- and three-degrees-of-freedom systems suggest a striking feature associated with these subtle 
bifurcations: the dimension typically exhibits a sharp, cusplike local minimum at the bifurcation. 
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The development of nonlinear dynamics has led to new un- 
derstandingof important physical processes, such as chaotic 
scattering [0, Q]. A scattering process is chaotic if it exhibits 
a sensitive dependence on initial conditions in the sense that, 
a small change in the input variables before the scattering can 
result in a large change in some output variables after the scat- 
tering. Apparently, in a scattering process, every physical tra- 
jectory is unbounded, but chaos arises because of a bounded 
chaotic set (chaotic saddle) in the scattering region where, for 
instance, the interaction between particles and potential oc- 
curs, as in a potential scattering system. Chaotic saddles are 
nonattracting and, dynamically, they lead to transient chaos 
Chaotic scattering is thus the physical manifestation of 
transient chaos in open Hamiltonian systems, which occurs in 
a variety of physical contexts [0] . 

An issue of fundamental importance in the study of chaotic 
scattering is to understand how the scattering dynamics 
changes characteristically as a system parameter is varied. 
There are three distinct classes of bifurcations in chaotic scat- 
tering: (1) routes to chaotic scattering from regular scatter- 
ing, (2) metamorphic bifurcations in which the chaotic saddle 
suddenly changes its characteristics, and (3) topological bifur- 
cations in which the topology of the chaotic saddle changes. 
Routes to chaotic scattering have been well documented, in- 
cluding abrupt bifurcations [0 |]] and saddle-center bifurca- 
tion followed by a cascade of period-doubling bifurcations 
There have also been efforts to study metamorphic bi- 
furcations in chaotic scattering, such as crisis [Q]. Topologi- 
cal bifurcations, characterized by topological changes in the 
subsets of infinite number of unstable periodic orbits (UPOs) 
embedded in the chaotic saddles, on the other hand, have not 
been well understood. One example of topological bifurca- 
tion is the massive bifurcation [||]. This type of bifurcations 
is subtle but it is expected to be important in high dimensions, 
as they can lead to distinct, physically observable scattering 
phenomena |||] . While there have been quantitative scaling re- 
sults concerning both routes to chaotic scattering [Q, [Io[ and 
metamorphic bifurcation [J7j], the understanding of topological 
bifurcations in chaotic scattering remains largely at the quali- 
tative level [B], particularly in high dimensions [0]. 

In this paper, we present a scaling law for the fractal di- 
mension of the chaotic saddle in a topological bifurcation that 
appears to be typical in both two- [pt] and three-degrees-of- 



freedom (DOF) Hamiltonian systems. We focus on the 
fractal dimension because of its importance in shaping phys- 
ically measurable quantities such as scattering functions, and 
because of the fact that most analytic scaling results associ- 
ated with various routes to chaotic scattering concern frac- 
tal dimension but, to our knowledge, there has been none in 
the context of topological bifurcations. Due to the subtlety of 
a topological bifurcation, obtaining quantitative scaling laws, 
even numerically, is highly nontrivial. To be concrete, we state 
our main result in the context of potential scattering, with the 
particle energy E as the bifurcation parameter. Let E c be the 
critical energy value at which a topological bifurcation occurs 
in the sense that, a class of infinite number of UPOs is de- 
stroyed and replaced by another. The principal result of this 
paper is then that, near the topological bifurcation, the fractal 
dimension of the set of singularities in any scattering function 
(to be described below) scales with the energy in the following 
way: 
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where D c is the dimension value at the bifurcation point. The 
remarkable observation is that the dimension-versus-energy 
curve exhibits a cusplike minimum at the bifurcation. We es- 
tablish Eq. (0) by a scaling argument, and we provide numer- 
ical verification using a class of 2-DOF scattering systems. 
We also argue that the scaling law should hold in 3-DOF sys- 
tems. Our result is interesting, as it represents the first, gen- 
eral, quantitative scaling result for topological bifurcation of 
chaotic scattering. 

To derive the scaling law (0), we first consider 2-DOF po- 
tential scattering systems. Our idealized, albeit representative, 
scattering system consists of three potential hills in the plane, 
as shown in Fig. The heights E\,E%, and E% of hills 1, 

2, and 3 satisfy E\,E2 > E3. Assume that hill 3 has a cir- 
cularly symmetric, quadratic maximum, and that it is close to 
the line connecting hills 1 and 2 so that the angle 9 in the trian- 
gle formed by the centers of all three hills is greater than 90°, 
as shown in Fig. 0(a). When the particle energy E is slightly 
above £3, a trapped particle can move back and forth between 
hills 1 and 2 following three different paths [H]: h, i, and j, 
as shown in Fig. 0(a). When the energy drops below E3, the 
path with greater deflection by crossing near the center of hill 

3, path j, is destroyed because of the appearance of a new for- 
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FIG. 1: A topological bifurcation in 2-DOF case: basic paths com- 
posing UPOs (a) before and (b) after the bifurcation. 
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FIG. 2: Construction of the Cantor set. 



bidden region around the center of hill 3, and is replaced by 
two new paths k, between hills 1 and 3 and between hills 2 
and 3, as shown in Fig. |j](b). For E w E 3 , the UPOs that 
form the skeleton of the chaotic saddle are thus composed of 
sequences of paths h, i, j, and k. A topological bifurcation [^] 
occurs at E c — E3, at which the class of an infinite number of 
UPOs containing path j is destroyed and replaced by another 
class of infinite number of UPOs that contain paths k. 

We focus on a fractal dimension that is physically most rel- 
evant: the dimension D of the set of singularities in a scat- 
tering function. Suppose that we launch particles toward the 
scatterer from a line outside the scattering region and measure 
a dynamical variable characterizing the outgoing particles af- 
ter the scattering {e.g., scattering angle), as a function of a 
variable characterizing the incoming particles before the scat- 
tering {e.g., impact parameter). Due to chaos, such a scatter- 
ing function typically contains an infinite number of singular- 
ities that constitute a fractal set, which is the set of intersect- 
ing points between the stable manifold of the chaotic saddle 
and the line [pi]]. Dynamically, the fractal set of singulari- 
ties is constructed by successive interactions of the particles 
with the potential hills. For example, consider the interval s 
of the line correspondingto particles that are first deflected by 
hill 1, as shown in Fig. 0. Three subintervals of this interval 
are then deflected by the other hills, with gaps between them 
corresponding to particles that are scattered to infinity right 
after the first interaction with hill 1. When E > E c , these 
three subintervals are defined by particles that go from hill 1 
to hill 2 along paths h, i, and j, respectively, as shown in Fig. 
^(a). When E < E c , while two of the subintervals are still 
defined by particles that follow paths h and i, the third one is 
now determined by particles that go from hill 1 to hill 3 and 
back to hill 1, following path k, as shown in Fig. ^(b). In the 
next step each one of the three surviving subintervals will split 
into three smaller subintervals, and so on. The splitting due to 
hill 3 is recorded at successive interactions with hills 1 and 2. 
This process defines a Cantor set with three linear scales cor- 
responding to the factors by which the intervals are reduced in 
each step: a w Sh/s, (3 w s;/s, and 7 « Sj/s (E > E c ) or 
7 sa Sfc/s {E < E c ), as shown in Fig. g. 

To obtain the dimension D of the Cantor set, we employ the 
following scaling argument. Suppose that the critical energy 



E c , the distance between the hills, and the effective radii of 
hills 1 and 2 at E c are of order unity. Near the bifurcation, the 
deflection angle <fr due to hill 3 of the particles corresponding 
to the edges of the intervals Sj and Sk, is smaller and greater 
than 90° for E > E c and E < E c , respectively. In both cases, 
I tan cj>\ is of order unity insofar as the angle 6 is close neither 
to 90° nor to 180° so that tan(7r - 6) ~ 0(1) [see Figs. 
and El]. Thus, setting | tan <j>\ to be of order unity provides an 
estimate of the width of the intervals Sj and Sfc, and hence an 
estimate of 7. For concreteness, say the maximum of the po- 
tential hill 3 (the bifurcation hill) has the following quadratic 
form: V = E c [l — (r/a) 2 ]. A straightforward calculation [ p"2| ] 
indicates that if | tan (f>\ ~ 0(1), then the impact parameter b 
scales with the energy difference: b s» b^\E — E c \ = boA, 
where bo is proportional to |tan^|. Therefore, the Cantor 
set is determined by the following scaling factors: a, f3, and 
7 70 A, where a and (3 are constants, and 70 is constant 
on each side of the bifurcation. Under the approximation 
that there is a self-similarity in the fractal structure, the di- 
mension D satisfies the following transcendental equation [Q]: 
a D + (3 d + 7( f A D = 1. For A and 5 = D — D(E C ) small, 
to the lowest order, this equation reduces to 

[A D iEc)] s /D{Ec) = c x sjom (2) 

where C is a positive constant. An asymptotic analysis of Eq. 
(Q) leads to [|l3]] our main scaling result ([!]) fll4]]. 

We now present numerical confirmation for the scaling law 
(jl|). To resolve the cusp structure implied by Eq. (|l]) is a 
highly nontrivial task. For numerical feasibility we thus con- 
sider a planar scattering system consisting of three nonover- 
lapping potential hills located at the vertices {xi,yi) (i — 
1,2,3) of an isosceles triangle [^[j. Each potential is repre- 
sented by the following function: Vt (x, y) — E.- L [1 — (r, /ai) 2 ] 
for r*j < ttj, where r, t = \/ (x — Xi) 2 + (y — yi) 2 , and the 
potential is zero elsewhere. The Hamiltonian is thus: H = 
|p| 2 /2m + V(x, y), where p is the two-dimensional momen- 
tum vector, and V(x,y) — 5Zi=i Vi{ x i V)- The advantage 
of utilizing quadratic potential hills is that the motion of the 
particle inside each potential can be solved exactly, making 
a high-precision dimension calculation possible. In our nu- 
merical experiment, we choose the following set of parameter 
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FIG. 4: A topological bifurcation in 3-DOF case: basic paths (a) 
before and (b) after the bifurcation. 



FIG. 3: Numerically obtained cusp scaling behavior in the fractal 
dimension about a topological bifurcation. 



values: E\ — E 2 — 10, E 3 = 1, ai = a 2 = 2, a 3 = 1, 
(xi,Vi) = -(x2,y 2 ) = (4,0), (x 3 ,y 3 ) = (0,2), and m = 1. 
A topological bifurcation thus occurs at E c = 1.0. We uti- 
lize the uncertainty algorithm [|l6|| to compute the dimension 
D of the set of singularities in scattering functions. For il- 
lustrative purpose, we choose the function to be the scatter- 
ing angle to infinity tp = tan -1 (p y /p x ) versus the loca- 
tion of the upward (initially) particle on the horizontal line 
at Vo = — 4. The result is shown in Fig. 0, in which the 
cusp behavior is unequivocal [ p"7[ ]. Linear least-squares fits 
between In [D(E) - D(E C )} and In \E - E c \ yield: D{E) - 
D(E C ) = (0.093 ± 0.005) x (1 - £)0-29±o.o3 for e < E c 
and D(E) ~ D(E C ) = (0.19 ± 0.01) x (E - l)0-3i±o.o3 for 
E > E c , as shown by the inset in Fig. [|. The scaling ex- 
ponents for both E > E c and E < E c agree well with the 
theoretically predicted one, the value of the dimension at E c , 
which we obtain numerically: D(E C ) = 0.314±0.005. In ad- 
dition, the proportional factor in Eq. (0) is greater for E > E c 
than for E < E c , which is also expected because, theoreti- 
cally, this factor is proportional to and 70 is greater in 
the former case. Overall, our computations lend strong cre- 
dence to the validity of the scaling ([j]). 

Our argument for the scaling law (0) can in fact be extended 
to chaotic scattering in 3-DOF Hamiltonian systems. As an 
example, we consider four potential hills but they are now lo- 
cated in the three-dimensional physical space at the vertices 
of a tetrahedron, as schematically illustrated by planarprojec- 
tion of the hills in Fig. Qja) for E > E c and in Fig. |(b) for 
E < E c . We focus on the dimension D of the intersecting 
set between the stable manifold of the chaotic saddle and a 
two-dimensional surface of the phase space, which is the di- 
mension of the set of singularities in a scattering function on 
two variables. Arguments similar to those in the 2-DOF case 
can then be used to derive the scaling law (0) near a topo- 
logical bifurcation, with the understanding that the dimension 



can in principle be either smaller or greater than 1 . When the 
radii of the potential hills are small, D(E C ) < 1 can arise, in 
which case dD(E) / dE diverges at E c and the dimension ex- 
hibits a cusp of the same nature as that in the 2-DOF case. If 
D(E C ) > 1, dD(E) I dE goes continuously to zero as E goes 
to E c , which implies that the dimension should be smooth. 

Can this smooth behavior be expected for a topological bi- 
furcation in 3-DOF scattering systems? To obtain an answer, 
we note the defining characteristic of a topological bifurca- 
tion: the topology of the UPOs is changed but from the stand- 
point of symbolic dynamics, no symbol is missing. In our 
3-DOF case it means that, before the bifurcation, all UPOs 
corresponding to periodic sequences of the nine basic paths 
shown in Fig. 0(a) must actually exist, while after the bifurca- 
tion, all those corresponding to sequences of the paths shown 
in Fig. ^(b) exist. In particular, the UPO represented by the se- 
quence {h, j} must exist when E > E c . Without loss of gen- 
erality, suppose that our 3-DOF system consists of three hard 
spheres of radius R at the vertices of a regular triangle with 
edge of unit length, and a smooth bifurcation hill located on 
the symmetrical axis perpendicular to the center of the trian- 
gle. In the inset of Fig. 0(a) is shown the periodic orbit {h,j} 
for E > E c in the limit E — > E c . A necessary condition for 
the existence of this orbit, and hence for the occurrence of the 
topological bifurcation, is then £ > 90° [see Fig. 0(a): the de- 
flection angle due to hill 3 is smaller than 90° when E > E c ], 
In the three-dimensional physical space, however, there is an 
additional geometric constraint: c > 1/2 V3, where c is the 
orthogonal distance from the center of the bifurcation hill to 
the lines connecting the hard spheres [Fig. 0(a)]. We obtain 
numerically that the conditions c > 1/2a/3 and ( > 90° are 
satisfied simultaneously only for R < 0.39. Numerical exper- 
iments indicate that, in this range of parameter values, D(E C ) 
is smaller than 1 . It is reasonable to assume that the same is 
valid when the hard spheres are replaced by smooth potential 
hills. This implies that for a generic topological bifurcation 
in 3-DOF scattering systems, the cusp behavior in Eq. (Q) is 
always expected 018]. 

In summary, our scaling analysis and numerical computa- 
tions reveal a striking behavior in the fractal dimension asso- 
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